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Abstract: Soil moisture is critical for vegetation growth in deserts. However, detailed data regarding the soil 
moisture distribution in space and time in the Gurbantungeut Desert of China have not yet been reported. 
In this study, we conducted a series of 7” situ observation experiments in a fixed sand dune at the southern 
edge of the Gurbantungegut Desert from February 2014 to October 2016, to explore the spatio-temporal 
variation of soil moisture content, investigate the impact of Haloxylon ammodendron (C. A. Mey.) Bungeon 
soil moisture content in its root zone, and examine the factors influencing the soil moisture spatial pattern. 
One-way analysis of variance, least significant difference tests and correlation analysis were used to analyze 
the data. The results revealed that the soil moisture content exhibited annual periodicity and the temporal 
variation of soil moisture content throughout a year could be divided into three periods, namely, a moisture- 
gaining period, a moisture-losing period and a moisture-stable period. According to the temporal and spatial 
variability, the 0—400 cm soil profile could be divided into two layers: an active layer with moderate variability 
and a stable layer with weak variability. The temporal variability was larger than the spatial variability in the 
active layer, and the mean profile soil moisture content at different slope positions displayed the trend of 
decreasing with increasing relative height and mainly followed the order of interdune area>west and east 
slopes>slope top. The mean profile soil moisture content in the root zone of dead H. ammodendron 
individuals was significantly higher than that in the root zones of adult and young individuals, while the soil 
moisture content in the root zone of adult individuals was slightly higher than that in the root zone of 
young individuals with no significant difference. The spatial pattern of soil moisture was attributable to the 
combined effects of snowfall, vegetation and soil texture, whereas the effects of rainfall and evaporation 
were not significant. The findings may offer a foundation for the management of sandy soil moisture and 
vegetation restoration in arid areas. 
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1 Introduction 


Soil moisture is the main factor restricting vegetation growth in arid and semi-arid areas (Hu et al., 
2009; Morenode et al., 2011), especially desert regions, and it not only affects plant individual 
development and growth types (Gad and Kelan, 2012) but also impacts plant species diversity, 
quantity and distribution (Kemp, 1983; Porporato et al., 2002). Soil moisture is also a critical 
component of water circulation in the soil—plant—atmosphere continuum (Cowan, 1965) and is the 
basis for the determination of vegetation carrying capacity and the management of sandy soil 
moisture (Rodriguez-Iturbe et al., 1999). However, soil moisture is highly variable in time and 
space across various scales, which we refer to as the spatio-temporal variability of soil moisture 
(Penna et al., 2009; Brocca et al., 2010; Heathman et al., 2012). This spatio-temporal variability 
exerts a strong influence on vegetation restoration (Meerveld and McDonnell, 2006; Chen et al., 
2007). Understanding this variability is crucial for analyzing the related hydrological processes 
(Hupet and Vanclooster, 2002; Choi and Jacobs, 2007) and landscape ecology (Penna et al., 2013), 
which is one of the major challenges in modern hydrology (Vereecken et al., 2008). 

Over the past several decades, the spatio-temporal patterns or variability of soil moisture have 
been extensively investigated by researchers in the fields of pedology and hydrology. These studies 
focused on various land use types, such as farmland (Baroni et al., 2013), forestland (Buttafuoco et 
al., 2005) and grassland (Hu et al., 2011), and different spatial scales, such as field scale (Heathman 
et al., 2012), watershed scale (Wilson et al., 2003; Heathman et al., 2012; Rosenbaum et al., 2012), 
continental scale (Li and Rodell, 2013) and global scale (Reichle et al., 2004). However, fewer 
studies have been conducted in arid environments, especially natural desert areas, owing to the 
difficulties and associated costs of sampling soil moisture (Gao et al., 2013; Zhang and Shao, 2013). 
Huang et al. (2015) analyzed the spatio-temporal distribution of soil moisture and salinity in the 
shelterbelt system of the Taklimakan Desert. Pan et al. (2015) compared the spatio-temporal 
variability of root-zone soil moisture in artificially revegetated and natural ecosystems in the 
Tengger Desert. Li et al. (2013) investigated the spatio-temporal variability of soil moisture at a 
desertified aeolian riparian ecotone on the Tibetan Plateau. 

The spatio-temporal variation of soil moisture is dependent on numerous factors, such as 
topography, soil properties, and vegetation. Jia and Shao (2013) reported that vegetation is the main 
factor influencing soil moisture patterns. However, Gomez-Plaza et al. (2001) suggested that 
although vegetation plays the major role in vegetated zones, soil texture and slope are the main 
regulators of soil moisture distribution in non-vegetated zones. Pan et al. (2008) also reported a 
difference in the influencing factors between dry periods and rainfall periods, where during rainfall 
events the main factors were local topography and vegetation, and during dry periods the soil 
texture was most important. Joshi and Mohanty (2010) indicated that soil texture and topography 
were two significant influencing factors affecting the spatio-temporal variation of soil moisture at 
point and remotely sensed footprint scales. However, the dominant factor is expected to vary under 
different hydrological settings (Western et al., 2004), and further studies of a variety of soils, 
terrains, and climatic conditions and over a wide range of scales are required (Zhang and Shao, 
2013). 

The Gurbantunggut Desert is the second largest desert and the largest fixed or semi-fixed desert 
in China (Qian and Wu, 2010). It is unique in terms of having been identified as a snow-dominated 
region among the four largest deserts in China, because the majority of its area lies north of 45°N, 
where the hydrological cycle is dominated by snowmelt (Barnett et al., 2005). There is a seasonal 
snow cover period that lasts for approximately 110 d (Zhou et al., 2012). However, the spatio- 
temporal variation of soil moisture content in fixed sand dunes has not been fully analyzed, and the 
rainfall, snowfall, evaporation, vegetation and soil texture may severely affect the spatio-temporal 
patterns of soil moisture. Therefore, it would be of great interest to examine the spatio-temporal 
variation of soil moisture and influencing factors in the Gurbantunggut Desert. 

In this work, a typical fixed sand dune on the southern edge of the Gurbantunggut Desert was 
selected as the study object, and in situ observations were performed from February 2014 to 
October 2016. The objectives of this study were to (1) evaluate the spatio-temporal variation of 
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sandy soil moisture, (2) assess the impact of Haloxylon ammodendron on the soil moisture content 
in the root zone, and (3) identify the factors influencing the spatial pattern of soil moisture. 


2 Materials and methods 


2.1 Study site 


The study site was located in the Beishawo experimental site (44°22'38"N, 87°55'13"E; 435 m 
a.s.l.) of the Fukang Desert Ecosystem Observation and Experiment Station, Chinese Academy of 
Sciences, on the southern edge of the Gurbantunggut Desert. The area of the study site has a 
temperate continental climate, with hot and dry summers in which the air temperature can reach 
44.2°C and cold winters in which the air temperature can drop to —42.2°C (Zhou et al., 2012). The 
annual mean temperature is 6.6°C, the mean annual precipitation is 70-180 mm, 25% of which is 
generally snowfall, and the annual pan evaporation is approximately 2000 mm (Xu et al., 2007). 
The snow cover period typically lasts 110-120 d from mid-November to mid-March of the 
following year, and the snow depth is approximately 20-30 cm. Ephemeral plants bloom and thrive 
in spring with the infiltration of snowmelt water and rainfall, when the soil moisture content is at 
its highest level (Dai et al., 2015). The interdune groundwater depth is 9-10 m. The soil is classified 
as aeolian soil (Arenosol, FAO), with a bulk density ranging from 1.54 to 1.65 g/cm?. The dominant 
plant species are H. ammodendron and Haloxylon persicum, and some annual species such as 
Ephedra distachya, Ceratocarpus arenarius and Artemisia desertorum also grow in the area. 


2.2 Experimental design and data collection 


Three soil moisture monitoring sections were established in a fixed dune on the southern edge of 
the Gurbantunggut Desert in December 2013. Two monitoring sections were set along the 
transverse section of the dune, and the third one, which was vertical to the above-mentioned two 
sections, was set in the interdune area. In total, 23 soil moisture monitoring points were set, and 
these were located on the west slope, slope top, east slope and interdune area (Fig. 1). The soil 
moisture monitoring points located in the interdune area were set in bare land and the root zone of 
H. ammodendron at various growth stages. According to the height, crown diameter and basal 
diameter data (Table 1), we classified the H. ammodendron plants as young, adult or dead 
individuals. Jn situ measurements of the soil moisture were performed from February 2014 to 
October 2016 in the soil profiles (0-400 cm soil layer) at 10-d intervals using neutron probes (Fig. 
1). At each monitoring point, three moisture content measurements were taken at depths of 10-390 
cm at 20-cm intervals. The mean was taken as the soil moisture content at each depth. The neutron 
probe was calibrated by comparison with data obtained via the oven-drying method. Specifically, 
we embedded an aluminum pipe in sand dune soils and then irrigated around the pipe. When the 
soil water infiltration reached a stable state, the scaling ratios (reading in soil/reading at soil 
surface) of the neutron probe at various depths were recorded. Soil samples were then taken from 
the same depths as the above-mentioned measurements and oven dried, and a correlational 
relationship between the soil moisture content and neutron probe scaling ratio was established. The 
calibration equations were as follows for 0—30 (Eq. 1) and 30—400 cm (Eq. 2) soil layers: 


0, =1.0916R,, R° =0.8398, (1) 
0, =0.6760R, — 0.0122, R° =0.9833, (2) 


where 6; is the soil volumetric moisture content (cm?/cm?) and Ro is the neutron probe scaling ratio. 

The influencing factors were measured at the west slope, slope top, east slope and interdune area 
in 2015. At each position, the rainfall was measured using four custom-built rain gauges. For each 
rainfall event, the mean of the values from the four rain gauges at each position was taken as the 
rainfall for that position. The evaporation was measured every other day using three custom-built 
micro-lysimeters at each position, and the mean value was taken as the evaporation for that position. 
In mid-March, five points at each position were selected to measure the snow depth, and the mean 
value of the snow depths at the five points was taken as the snow depth for that position. The 
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vegetation coverage and density were measured in June using the traditional quadrat survey 
method. The soils were sampled in 0-200 cm soil layers for each position and the particle 
composition was measured using a laser particle size analyzer (Mastersizer 2000, Malvern 
Panalytical, UK). 


Table 1 Morphological features of Haloxylon ammodendron in different growth stages 


Growth stage Height (m) Crown diameter (mxm) Basal diameter (cm) 
Young (n=3) 1.7-2.1 1.5x1.3-1.7x1.5 5.4-7.3 
Adult (n=4) 3.6—4.2 1.9x3.0—4.2x4.0 11.5-17.2 
Dead (n=1) 2:2 2.4x2.3 8.0 


Note: n, the number of plants. 
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Fig. 1 Distribution of soil moisture observation sites in the west slope, slope top, east slope and interdune area. 
(a) and (b) are the two sections along the direction from West to East, and (c) is the bird view of the horizontal 
distribution of observation sites. Capital letters and numbers are the code of observation sites. 


2.3 Data analysis 


The coefficient of variation (CV) was selected to characterize the variability of soil moisture, and 
the CV values were analyzed according to the method suggested by Nielsen and Bouma (1985). 
The magnitude of the variability was considered weak when CV<10%, moderate when 
10%<CV<100% and strong when CV>100%. The spatial CV was calculated from the mean value 
and standard deviation for different monitoring points at each position. For the temporal CV, we 
first calculated the mean soil moisture content at each position at one time point, then calculated 
the mean value and standard deviation of all time points in a year. We used the mean profile soil 
moisture content to assess the differences in soil moisture content at different positions and for 
different growth stages of plants, which was obtained by calculating the mean soil moisture content 
at each depth during a year and then calculating the mean soil moisture content in a profile (0-400 
cm soil layers). For the influencing factors, the cumulative soil evaporation capacity was the sum 
of each measurement of soil evaporation, and the cumulative rainfall capacity was defined 
similarly. Differences in the soil moisture content at different positions and in the root zone of H. 
ammodendron at different growth stages were determined using SPSS 16.0 (SPSS Inc., USA). 
Graphs were plotted using Origin 8.0 (OriginLab Corporation, USA). 
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3 Results and discussion 


3.1 Temporal variation of soil moisture content 


On the basis of the soil moisture measurements during the three study years, the temporal variation 
of the soil moisture content throughout the year could be divided into three periods, namely, a 
moisture-gaining period (from March to mid-May), a moisture-losing period (from mid-May to 
October), and a moisture-stable period (from November to February of the following year). Figure 
2 shows the temporal variation of soil moisture content during 2014-2016, revealing that the soil 
moisture content exhibited a similar temporal variation pattern each year. As we can see, in early 
spring the soil moisture content was low and almost invariant throughout nearly the entire soil 
profile, except at the soil surface. Meanwhile, the soil moisture content near the soil surface 
remained consistently high on account of the snow cover and some degree of freeze—thaw action. 
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Fig. 2 Temporal variations of soil moisture content in 2014 (a), 2015 (b) and 2016 (c). Values are means of all 
observation sites (n=21). 


From March, the soil moisture content in the shallow layer (0-50 cm) started to increase 
significantly with snowmelts and the rising temperature and remained at a high level during the 
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whole of spring (from March to May); however, the soil moisture content in the deep layer (S0— 
400 cm) changed very little. In June, the soil moisture content began to fall rapidly in the surface 
soil layer owing to the intense solar radiation and the formation of a dry layer on the sand dune 
surface. The soil moisture content remained continuously low during the entire summer, except 
during rainfall events, which resulted in temporary increases. Meanwhile, the wetting front moved 
down in a step-by-step manner in the deep soil layer and the soil moisture content decreased 
gradually during the entire summer, although the soil moisture content was still relatively higher in 
the summer and autumn seasons. During winter, the soil moisture content decreased to a low level 
and remained almost invariant. 

These results are in agreement with those reported by Huang et al. (2015), who found that the 
soil moisture content (0—10 cm) in the Tarim Desert Highway Shelterbelt was the highest in spring 
and then decreased over time owing to evaporation and considered temperature to be the main 
influencing factor. Similarly, Li et al. (2010) compared the difference in soil moisture (0-120 cm) 
temporal variation patterns between shrubs and bare land in the Gurbantunggut Desert and found 
that the topsoil moisture content was the highest in April and then decreased over time until 
September, while the moisture content in the upper soil layers was higher than that in the deeper 
soil layers during spring and the opposite trend was observed in summer and autumn. These authors 
attributed the increased soil moisture content in the upper layers in spring to the infiltration of 
melting water and the decreased soil moisture content in the upper layers in summer to evaporation. 
In our study, we also discovered a similar temporal variation of soil moisture; however, there were 
some differences in the soil moisture content and soil moisture change time point owing to the 
differences in climate, vegetation, underlying surface and so on. 

In the Gurbantunggut Desert, the soil moisture content is influenced by snowmelt water 
infiltration recharge, precipitation infiltration recharge, plant absorption, and evapotranspiration. 
Snowmelt water is an important water source that accounts for 18.7%-—30.0% of the annual 
precipitation (Li, 1991; Xu et al., 2007; Zhou et al., 2012). The soil moisture content increases 
rapidly in spring owing to melting of the snow on the sand dune surface (Zhou et al., 2012). In mid- 
May, the temperature rises rapidly and soil evaporation increases, while shrubs begin to bud and 
large numbers of ephemeral plants grow, causing a net loss of soil moisture. By June, H. 
ammodendron has reached a vigorous growth stage and a large number of herbaceous plants grow 
vigorously, which, in addition to the high temperature, cause the consumption of soil moisture to 
increase greatly. From May to October, the soil moisture content shows a downward trend, despite 
brief increases due to rainfall events. In early November, the sand dune surface is already covered 
by snow, the temperature is very low and the moisture in the soil has already frozen, while the H. 
ammodendron plants have also disrobed of their leaves and consumed minimal soil moisture, and 
consequently the soil moisture content changes little during the entire winter. 


3.2 Variation of soil moisture content with slope position 


Figure 3 presents the variation of the mean profile soil moisture content with slope position. In 
2014, the difference in soil moisture content was not significant between the west slope and east 
slope, or between the east slope and slope top, although it was significant between the west slope 
and slope top. The differences in soil moisture content between the interdune area and the above- 
mentioned three positions were significant, and the soil moisture content in the interdune area was 
relatively higher. Specifically, the soil moisture content in the interdune area was 1.08, 1.15 and 
1.10 times that in the west slope, slope top and east slope, respectively. In 2015, the difference in 
soil moisture content was not significant between the slope positions, i.e., the west slope, slope top 
and east slope, but it was significant between the interdune area and the east slope and slope top. 
The difference in soil moisture content between the interdune area and west slope was not 
significant. However, the soil moisture content in the interdune area was still higher than that in the 
other slopes. Specifically, the soil moisture content in the interdune area was 1.09, 1.14 and 1.09 
times that in the west slope, slope top and east slope, respectively. In 2016, the soil moisture content 
in the slope top was the lowest and displayed a significant difference compared with the other 
positions. The difference was not significant among the west slope, interdune area and east slope; 
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however, the soil moisture content in the interdune area was still the highest. Specifically, the soil 
moisture content in the interdune area was 1.04, 1.18 and 1.01 times that in the west slope, slope 
top and east slope, respectively. The three-year-average soil moisture content was highest in the 
interdune area, and the soil moisture content in this area was significantly different from that in the 
slope top but not significantly different from that in the two slopes. The difference in the soil 
moisture content was not significant between the west slope and east slope. Meanwhile, the soil 
moisture content in the slope top was the lowest and also not significantly different from that in the 
two slopes. Overall, the soil moisture content at different slope positions displayed the trend of 
decreasing with increasing relative height, and mainly followed the order of interdune area>west 
and east slopes>slope top. 
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Fig. 3 Variances of mean profile soil moisture content in west slope, slope top, east slope, and interdune area in 
2014 (a), 2015 (b), 2016 (c) and the average for the period 2014—2016. Bars indicate standard errors. West slope, 
east slope and interdune area, n=4; slope top, n=2. Different lowercase letters denote significant difference among 
different positions at P<0.05 level. 


Other researchers have also studied the spatial patterns of soil moisture in deserts. For example, 
Pan et al. (2009) investigated the spatial pattern of soil moisture and influencing factors in the 
Tengger Desert and found that the soil moisture content was negatively correlated with relative 
elevation. Zhou et al. (2013) also demonstrated that the soil moisture content in the slope toe was 
the highest in the Gurbantunggut Desert. These results were all consistent with our finding that the 
soil moisture content decreased with increasing relative elevation. 

On account of plants being relatively dense in the interdune area, the degree of cover was greater, 
and the "wet island" effect (Yang et al., 2011) was obvious. Taking the soil particle size composition 
into consideration, Pan et al. (2009) reported the existence of a positive correlation between the 
soil moisture content and clay or silt content. The soil fine particle (<0.100 mm) content in the 
interdune area was relatively higher than that at other positions. Furthermore, snowmelt water flows 
along the slope surface to the interdune area in spring to afford an increased soil moisture content. 
All of these reasons account for the relatively higher soil moisture content in the interdune area. 
Few plants grow on the slope top, which consists mostly of bare sand, and has a limited capacity 
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for retaining water, so that the snowmelt water flows away during snowmelt periods. Consequently, 
the soil moisture content was the lowest at the slope top. 


3.3 Temporal and spatial variability of soil moisture content with soil depth 


Figure 4 shows the variation of the mean annual soil moisture content and temporal and spatial 
coefficient of variations with soil depth. The soil moisture content exhibited an obvious vertical 
variation and was higher and decreased rapidly with increasing depth in the 0—50 cm soil layer. At 
the same time, the soil moisture content was relatively constant with depth in the 50—400 cm soil 
layer. The temporal variability of the soil moisture content was evaluated using the temporal CV, 
which decreased with increasing soil depth in the 0-100 cm soil layer (Figs. 4b, e and h). The 
temporal variability was assigned as moderate in this layer according to Nielsen and Bouma (1985). 
In the 100—400 cm soil layer, the temporal CV was small and relatively stable and the temporal 
variability was weak. Similarly, the spatial CV decreased with increasing soil depth in the 0-50 cm 
soil layer with moderate variability (Figs. 4c, f and i) and the spatial CV was small and relatively 
stable in the 50-400 cm soil layer with weak variability. According to the temporal and spatial 
variability, the 0-400 cm soil profile could be divided into two layers: an active layer with moderate 
variability and a stable layer with weak variability. Specifically, the soil moisture content exhibited 
moderate variability in the 0-100 cm layer in terms of the temporal CV and in the 0-50 cm layer 
in terms of the spatial CV. Furthermore, the temporal variability was greater than the spatial 
variability in the active layer. 
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Fig. 4 Variations of mean annual soil moisture content, and temporal and spatial coefficients of variation with 
soil depth in 2014 (a, b, c), 2015 (d, e, f) and 2016 (g, h, i). West slope, east slope and interdune area, n=4; slope 
top, n=2. 


The temporal variability of soil moisture with soil depth was also reported by Pan et al. (2015), 
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who compared the differences in soil moisture in the root zone between artificial and natural 
vegetation in the Tengger Desert and found that the temporal variation was more obvious in 0—60 
cm soil layers than in 0-300 cm soil layers, while the soil moisture in 0—60 cm soil layers was more 
easily changed than that in 0-300 cm soil layers. Furthermore, Gao et al. (2015) investigated the 
spatio-temporal variability of soil moisture in the Loess Plateau and discovered that the temporal 
CV decreased with increasing soil depth in 0-160 cm soil layers and was essentially constant in 
160-300 cm soil layers. These results were all similar to our findings, despite some differences in 
the soil profile division. 

With respect to the spatial variability of soil moisture with depth, Wang et al. (2013) explored 
the vertical distribution of soil moisture in the Loess Plateau and found that the spatial CV 
decreased with increasing soil depth in 0—60 cm soil layers, fluctuated in 60-200 cm soil layers, 
and was constant in 200—400 cm soil layers. This is consistent with our results for upper soil layers, 
but inconsistent with our results for deeper soil layers in terms of the fluctuation of the CV. This 
discrepancy was mainly attributable to the differences in spatial scale, i.e., the microtopography 
scale in our study and the regional scale of Wang et al. (2013). A larger scale may be associated 
with more variability (Western and Bléschl, 1999), such as heterogeneity in soil texture, vegetation, 
terrain and so on, all of which can affect the soil moisture distribution. 

Soil moisture content is significantly affected by precipitation, soil surface evaporation, wind, 
heat conduction and soil pore water diffusion (Famiglietti et al., 1998; Li et al., 2013). In our study 
area, as the temperature rose in spring, snowmelt water infiltrated into the sand dune soil, and there 
was a frozen soil layer under the sand dune surface that impeded snowmelt water infiltration (Li et 
al., 2010). These phenomena jointly account for the high mean annual soil moisture content in the 
shallow soil layer. In terms of temporal variability, the active layer was more easily affected by 
meteorological factors, vegetation growth and so on. For example, rainfall resulted in a sharp 
increase in soil moisture content and the soil moisture content dropped rapidly with evaporation, 
while the soil moisture content in the active layer changed greatly owing to snowmelt water 
infiltration in spring. Vegetation growth also contributed to the difference in soil moisture content 
in different seasons. The spatial variability mainly resulted from the spatial variability of 
vegetation, soil texture and terrain. However, the spatial variability of these influencing factors was 
less than the temporal variability in the arid zone, such that the temporal variability of soil moisture 
was greater than the spatial variability in the active layer. In the stable layer, these factors exerted 
less of an influence on the soil moisture, and consequently the soil moisture variability was weak 
and the CV remained relatively constant. 


3.4 Variation of soil moisture content in the root zone of H. ammodendron at different 
growth stages 


Figure 5 presents the variation of the mean profile soil moisture content in the root zone of H. 
ammodendron at different growth stages. In 2014 and 2016, the difference in the mean profile soil 
moisture content in the root zone of plants was not significant between young and adult individuals, 
although it was significant between dead individuals and the other two growth stages. Similar 
results were observed in 2015, except that the difference between adult and dead individuals was 
not significant. The mean profile soil moisture contents in the root zones of young, adult, and dead 
individuals were 0.035, 0.037 and 0.041 cm3/cm?, respectively, in 2014; 0.033, 0.036 and 0.039 
cm?/cm?, respectively, in 2015; and 0.034, 0.038 and 0.043 cm?/cm’, respectively, in 2016. The 
three-year average mean profile soil moisture content in the root zone of dead individuals was 
significantly higher than that for adult and young individuals, while the soil moisture content in the 
root zone of adult individuals was slightly higher than that for young individuals with no significant 
difference. 

A number of studies have demonstrated that the soil moisture content in shrubland is higher than 
that in bare land, especially in arid and semi-arid environments (D'Odorico et al., 2007), owing to 
a combination of the "canopy effect" and "soil effect" (Gomez-Aparicio et al., 2005). In our study 
area, the soil moisture content in the root zone of H. ammodendron was affected by stemflow water 
collection (Wang et al., 2013), evapotranspiration (Hu et al., 2011) and tree canopy shading. For 
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the dead individuals, there was only stemflow water collection and leaf transpiration did not occur; 
at the same time, cracks formed between the trunk and soil owing to withering and cracking of the 
trunk, such that rainwater could easily infiltrate into the root zone soil to afford a higher soil 
moisture content. In terms of the adult and young individuals, the adult individuals possess a larger 
canopy that allows the collection of more rainwater, as reported by Yuan et al. (2016), who found 
that plants with more small branches and higher leaf biomass were more contribute to the formation 
of stemflow. Meanwhile, the tree canopy shading effect of adult individuals is more efficient. With 
respect to evapotranspiration, even if adult individuals lose more water via transpiration, this effect 
may be alleviated by canopy shading. Consequently, these effects jointly account for the higher soil 
moisture content in the root zone of adult individuals than in that of young individuals. 
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Fig. 5 Variances of mean profile soil moisture content in the root zone of Haloxylon ammodendron at different 
growth stages in 2014 (a), 2015 (b), 2016 (c), and the average of the period 2014-2016 (d). Bars indicate standard 
errors. Dead individual, n=1; adult individual, n=4; young individual, n=3. Different lowercase letters denote 
significant difference among different growth stages (P<0.05). 


3.5 Factors influencing spatial patterns of soil moisture content 
3.5.1 Rainfall 


Rainfall has a low spatial variability on a small regional scale, and the influence of rainfall on soil 
moisture content spatial patterns mainly depends on whether the amount of rainfall is sufficiently 
high to form surface flow or soil lateral seepage. Figure 6 shows the overall distribution of rainfall 
events and the variation of mean rainfall (individual rainfall event) with slope position. During the 
experimental period (from April to October) in 2015, 19 rainfall events were observed. The total 
rainfall was approximately 110.0 mm and most fell between April and September. The cumulative 
rainfall on the west slope, slope top, east slope and interdune area was 104.7, 105.5, 114.7 and 
111.4 mm, respectively, without significant differences. In terms of single rainfall events, the largest 
rainfall was 14.4 mm and the mean rainfall was 5.7 mm. The average soil infiltration rate in bare 
sandy land was 472.2 mm/h, which is far greater than the maximum rainfall intensity of 14 mm/h 
(Liu et al., 2015), such that surface runoff did not occur during natural rainfall for either the sand 
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dune or interdune area (Dai et al., 2014) and the rainfall water mainly infiltrated vertically. 
Consequently, rainfall did not exert a significant influence on the soil moisture content spatial 
pattern. 
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Fig. 6 Temporal variation of rainfall events (a) and variances of mean rainfall (individual rainfall) in west slope, 
slope top, east slope and interdune area (b). Bars indicate standard errors (n=4). 


3.5.2 Soil evaporation 

Figure 7 shows the variation of the cumulative soil evaporation (from May to September) with 
slope position. There was a clear linear correlation between the monthly evaporation capacity and 
monthly rainfall. Soil evaporation was mainly affected by rainfall. The average cumulative soil 
evaporation capacity was approximately 143.9 mm, and the cumulative soil evaporation capacities 
for the west slope, slope top, east slope and interdune area were 143.6, 142.3, 144.7 and 145.2 mm, 
respectively. The variance analysis results revealed no significant differences in the cumulative soil 
evaporation capacity with the slope position. Therefore, soil evaporation did not exert a significant 
influence on the soil moisture content spatial pattern. 
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Fig. 7 Variances of cumulative soil evaporation capacity (from May to September) in west slope, slope top, east 
slope and interdune area. Bars indicate standard errors (n=3). The error bar of soil moisture in west slope is too 
short and thus not shown. 


3.5.3 Snowfall 
Figure 8 shows the variation of the snow depth at the various slope positions. No significant 
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difference in the snow depth was observed between the west slope and slope top or between the 
east slope and interdune area. However, the snow depths at the east slope and interdune area were 
significantly greater than those at the west slope and slope top. As the prevailing wind direction at 
the study site in winter is northwest, the west slope was the windward slope and the east slope was 
the leeward slope, such that the snow on the west slope and slope top was blown away and piled 
up on the east slope and in the interdune area, causing the snow depths at the east slope and 
interdune area to be greater than those at the west slope and slope top. Furthermore, during the 
snow melting period, snowmelt water infiltration affected the spatial distribution of sand dune soil 
moisture content. As the snowmelt process was faster than the frozen soil layer melting process, 
the deep frozen soil layer impeded the snowmelt water infiltration until total melting of the frozen 
soil layer had occurred (Li et al., 2010). The snowmelt water therefore flowed along the slope 
surface to the interdune area. Therefore, snowfall exerted a significant influence on the soil 
moisture content spatial pattern. 
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Fig. 8 Variances of snow depth in west slope, slope top, east slope and interdune area. Bars indicate standard 
errors. Different lowercase letters denote significant difference (P<0.05) among different slope positions. 


3.5.4 Vegetation 

Vegetation influences soil moisture variability in several ways, such as the canopy changing the 
throughfall and shading the land surface to affect the rate of evaporation, the root activity and leaf 
litter affecting soil hydraulic conductivity, and the roots extracting moisture for transpiration from 
the soil profile (Famiglietti et al., 1998; Pan et al., 2008). As H. ammodendron is the dominant 
species in the Gurbantunggut Desert, we measured the coverage and density of H. ammodendron 
with slope position (Fig. 9). The vegetation coverage and density in the interdune area were higher 
than those at the west slope, slope top and east slope. Specifically, the vegetation coverage in the 
interdune area was 1.29, 1.20 and 1.14 times that at the west slope, slope top, and east slope, 
respectively, while the vegetation density in the interdune area was 1.70, 1.89 and 1.61 times that 
at the west slope, slope top and east slope, respectively. The correlation analysis results revealed 
that the soil moisture content was positively correlated with the vegetation density (P<0.05), as 
also reported by Pan et al. (2009), who found that the soil moisture content was positively correlated 
with shrub coverage. In semi-arid and arid environments, especially in natural desert areas, 
vegetation has a substantial impact on rainfall redistribution, and 5%-10% of rainfall is directly 
transported to the root zone soil by stemflow (Navar and Bryan, 1990; Mauchamp and Janeau, 
1993; Martinez-Meza and Whitford, 1996). The root channels and soil pores allow water to more 
easily infiltrate deep soil layers, such that plants can use the soil moisture in dry periods and soil 
moisture evaporation decreases (Nulsen et al., 1986; Reynolds et al., 1999). At the same time, the 
"shading effect" of H. ammodendron can reduce soil evaporation. Chen et al. (2014) investigated 
the difference in soil evaporation in the root zone of H. ammodendron and found that soil 
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evaporation increased gradually with increasing distance to the H. ammodendron base stem. This 
increased soil evaporation further reduced the soil moisture content, exacerbating the difference in 
soil moisture content distribution, as also reported by Li et al. (2010). Therefore, vegetation exerted 
a significant influence on the spatial pattern of soil moisture content. 
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Fig.9 Variances of vegetation coverage and density in west slope, slope top, east slope and interdune area 


3.5.5 Soil particle size composition 


Soil texture influences the hydrophysical properties and movement of soil moisture and therefore 
its distribution characteristics. Table 2 shows the variation of the soil particle size composition with 
slope position. The interdune area contained significantly greater proportions of soil particle sizes 
of 0.000—0.002, 0.002—0.050 and 0.050—0.100 mm compared with the west slope, slope top and 
east slope, whereas the slope top contained a significantly smaller proportion of soil particles of 
size 0.050—0.100 mm compared with the west and east slopes. The interdune area also contained 
significantly lower proportions of soil particle sizes of 0.100—0.250 and 0.250—0.500 mm than the 
west slope, slope top, and east slope, while the slope top contained significantly greater proportions 
of soil particles of size 0.250—0.500 mm. The proportions of soil particles of size 0.500—1.000 mm 
was very low at all slope positions without significant differences. Overall, the soil in the interdune 
area contained more fine soil particles and the soil at the slope top contained more large soil 
particles. We also calculated correlation coefficients between the mean soil particle size ranges and 
soil moisture content and found significant linear positive correlations between the soil particle 
sizes of <0.002, 0.002—0.050, and 0.050—0.100 mm and the soil moisture content and a significant 
linear negative correlation between the soil particle size of 0.100—0.250 mm and the soil moisture 
content; no significant correlation was observed for the soil particle sizes of 0.250—0.500 mm and 


Table 2 Soil particle size composition in west slope, slope top, east slope and interdune area and correlation 
coefficient between soil particle size and soil moisture content 


Soil particle size composition (%) 
Site 


<0.002 mm 0.002—0.050 mm 0.050—-0.100 mm 0.100—-0.250 mm 0.250—0.500 mm 0.500—1.000 mm 


West slope 0.063° 3.185° 20.424° 65.233* 11.094° 0.000° 
Slope top 0.023> 2.978" 9.393° 66.512* 21.077° 0.018" 

East slope 0.166" 3.404° 18.284” 62.663? 15.472° 0.010% 
Interdune area 0.690* 5.481° 29:532% 57.327" 6.9694 0.002° 
Correlation coefficient 0.956" 0.958" 0.975" -0.963" —0.917 —0.703 


Note: Data in the table is the mean value of soil particle size content in 0-200 cm soil layers. Values within a column followed by different 
lowercase letters are significantly different among different slope positions at P<0.05 level. * denotes significant correlation between soil 
particle size and soil moisture content at P<0.05 level. 
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0.500—1.000 mm. Similar findings were reported by Pan et al. (2008), namely, positive correlations 
between the soil moisture content and soil clay and silt contents and a negative correlation between 
the soil moisture content and sand content, as also found by Gao et al. (2011). Therefore, the soil 
particle size composition exerted a certain influence on the soil moisture content spatial pattern. 


4 Conclusions 


In this study, we investigated the temporal variation and spatial distribution of soil moisture and 
analyzed its spatio-temporal variability and the corresponding influencing factors. We found that 
the soil moisture content exhibited annual periodicity, where the temporal variation of soil moisture 
content throughout a year could be divided into three periods, namely, a moisture-gaining period, 
a moisture-losing period and a moisture-stable period. According to the temporal and spatial 
variability, the 0-400 cm soil profile could be divided into two layers: an active layer with moderate 
variability and a stable layer with weak variability. Meanwhile, the soil moisture temporal 
variability was larger than the spatial variability in the active layer. The mean profile soil moisture 
content at different slope positions displayed the trend of decreasing with increasing relative height 
and mainly followed the order of interdune area>west and east slopes>slope top. The mean profile 
soil moisture content in the root zone of dead H. ammodendron plants was significantly higher than 
that in the root zones of adult and young individuals, while the soil moisture content in the root 
zone of adult individuals was slightly higher than that in the root zone of young individuals with 
no significant difference. The spatial pattern of soil moisture was attributable to the combined 
effects of snowfall, vegetation, and soil texture, whereas the effects of rainfall and evaporation were 
not significant. Our study may offer a foundation for the management of sandy soil moisture and 
vegetation restoration in arid areas and provide a valuable reference for the fields of soil and water 
conservation, desertification control, and arid zone hydrology. 
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